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Abstract 

We present a model-independent analysis of K + —> 7 t + £ + £~ and Ks —> n°£ + £~ 
decays, including K —>• 3tt unitarity corrections and a general decomposition of the 
dispersive amplitude. From the existing data on K + —> 7r + e + e~ we predict the 
ratio R = B(K + —*■ tt + /i + p - )/ B(K + —> n + e + e~) to be larger than 0.23, in slight 
disagreement with the recent measurement R = 0.167 ± 0.036. Consequences for 
the —> 7r ± e + e _ charge asymmetries and for the Kj j —>■ 7r°e + e _ mode are also 

discussed. 


* Work supported in part by TMR, EC-Contract No. ERBFMRX-CT980169 
(EURODA<PNE) 


1 Introduction 


Radiative nonleptonic kaon decays continue to provide interesting information on 
the structure of the weak interactions at low energies. Among them, the flavour¬ 
changing neutral current (FCNC) transitions K —> n£ + £~, induced at the one- 
loop level in the Standard Model, are well suited to explore its quantum structure 
and, possibly, its extensions a a- On the experimental side, the high-precision 
measurements already achieved at Brookhaven (AGS) |3j and Fermilab (KTeV) [Q, 
and foreseen in the near future at Frascati (KLOE) ||, call for a thorough theoretical 
investigation. 

The K + —> 7 t + £ + £~ and Kg —> n°£ + £~ (£ = e,/i) channels, expected to be 
dominated by long-distance dynamics through one-photon exchange (K —* 717 *), 
were studied at leading order in the chiral expansion in Ref. [||J. Once the unknown 
combination of local terms arising in this framework is fixed by comparison with 
K + —> n + e + e~ data ||, the process K + —> 7 r + /r + /i~ can be predicted without further 
assumption^. Recently, this last mode has been observed [0], yielding a value for 
R = r(K+ —> 7 r + /i + /i _ )/r(A' + —> 7 r + e + e~) that is ~ 2 a below the leading-order 
prediction. This experimental result has motivated the present study of K —> tt'j* 
form factors beyond the leading order in the chiral expansion. We have carried out 
this program by including unitarity corrections from K —> mm together with the 
most general polynomial structure consistent with the chiral expansion up to 0(p 6 ). 
In this way we perform a model-independent analysis of the existing K + —> n + £ + £~ 
data leading to a consistent fit of both rate and dilepton invariant-mass spectrum 
of the electron channel. On the other hand, we demonstrate that the persisting 
discrepancy in the ratio R cannot be accommodated in the Standard Model. 

Besides this phenomenological analysis, we present a general discussion of theo¬ 
retical predictions for the polynomial part of the K + (Kg) —> 7T + (7r°)7* form factors. 
We show that it is very difficult at present to estimate the branching ratios of 
the Kg —> 7r°£ + £ - modes without strong model-dependent assumptions. More ac¬ 
cessible are the dilepton invariant-mass spectrum and, correspondingly, the ratio 
T(K S —> n°n + /i~)/T(Kg —> 7T°e + e~). 

We conclude our analysis by discussing the impact of the long-distance K —» 
7 T7* —» 7re + e“ transitions for estimating CP -violating observables. In particular, 
we analyse the possibilities for disentangling direct and indirect C'P-violating com¬ 
ponents of the Kg —> 7r°e + e _ amplitude. Moreover, we estimate the effect of the 
unitarity corrections on the charge asymmetry of Ii ± —> 7r ± e + e _ decays. 

Whe neutral channels (still unmeasured) require additional model-dependent assumptions as 
discussed in Sect. |j. 
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The plan of the paper is the following: in the next section we present the model- 
independent analysis of K —» 7ry* —> tt£ + £~ transitions. Section |3] is devoted to 
explore the physics behind the polynomial coefficients of the K —> 7ry* form factors. 
We turn to a discussion of CP violating observables in Sect, ft]. Our main conclusions 
are summarized in the last section. 


2 Model—independent analysis 

2.1 The FCNC transitions K —>■ i t£ + £~ {£ = e,/i) are dominated by single virtual 
photon exchange (K —> nj* —> tt£ + £~) if allowed by CP symmetry. This is the case 
for K + and Kg decays where the amplitude is determined by an electromagnetic 
transition form factor in the presence of the nonleptonic weak interactions: 


i 


d 4 xe^(n(p)\T{J^(x)C AS=l (0)}\K(k)) 


W(z ) 


z(k +pY - (1 - rl)q M 


k 2 = Mk , p 2 = M 2 , q = k-p , z = g 2 / M l > W = M n /M K , 


( 1 ) 


where £as=i is the strangeness changing nonleptonic weak Lagrangian and J^ m is 
the electromagnetic current. The dynamics of the decays is completely specified by 
the invariant functions W + (z) and Ws(z). As dictated by gauge invariance, these 
functions vanish to lowest order in the low-energy expansion |I]. To account for 
this chiral suppression we have pulled out a factor 1/(47r) 2 in the definition ([[]). 
Since the Cabibbo angle is approximately compensated by the nonleptonic octet 
enhancement, the natural magnitude of W(z) is expected to be GpM\. 

With these conventions^] the decay amplitude takes the form 


A(n(k) ^ n(p)e + ( p+ )r(pz) 


M 2 (a v w (z)(k + pYuiipJ)^Vtip + ) . (2) 

m k\ 


The spectrum in the dilepton invariant mass is then given by 


dT 

dz 


a 2 M, 


K 


127t(47t) 


A 3/2 (l ,z,rl) 


4^ I 1 + 2^-J \W(z)\- 


( 3 ) 


with re = me/M K and 4r 2 < z < (1 — r 7T ) 2 . 

2 In Ref. ^ the form factor V(z) = W(z)/G%M% was used instead of W(z). Here we prefer 
to avoid introducing the coupling constant G% in the model-independent analysis. We use the 
occasion to point out two misprints in Eq. (4.27) of Q; the correct formula is given by (jdj). 
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Figure 1: K —> 3n contribution to the effective K —> nj* vertex. 

2.2 The form factors Wi(z) {i = +, S) are analytic functions in the complex z— 
plane cut along the positive real axis. The cut starts at z — 4 r 2 with the two-pion 
threshold. For the small dilepton masses occurring in the decays, one expects the 
7 r + 7 r _ intermediate state to play the dominant role in the dispersion relations for 
Wi(z). The contribution of higher-mass intermediate states can be described by a 
low-order polynomial in z. This is known to be a very good approximation for the 
K + K~ contribution, for instance |IJ. 

We therefore decompose the form factors IF 7 ,; {z) as 

Wi{z) = G F M 2 K Wf°\z) + Wr(z) , (4) 

where W™(z) denotes the contribution from the two-pion intermediate state. To 
leading nontrivial order in the chiral expansion |l|, the W- >ol (z) are constants (but 
for a negligible contribution due to the kaon loop). The functions W^^z) can be 
calculated from the diagram shown in Fig. [Tj. where the K —> 37 T vertex is taken 
from the leading nonleptonic weak Lagrangian of Oijp 2 ). 

Here we go one step further. We use the physical K —>• 3tt amplitude expanded 
up to 0(p 4 ) (as required by present K —► 37 r data) and include the electromagnetic 
form factor F(z) (normalized to F( 0 ) = 1 ) for 7 r + 7 r _ —> 7 * to first nontrivial order. 
The relevant K —> 37 t amplitudes are expanded as [|[ 

A{K + ^ tt+tt+tt-) = 2a c + (b c + b 2 )Y + 2c c (y 2 + X 2 /3) + (d c + d 2 )(Y 2 - X 2 /3 ) , 
A{Ks ^ TT+TT-TT 0 ) = h 2 X-p 2 XY, (5) 

with 

Si = (k- Pi ) 2 , s 0 = i(si + s 2 + s 3 ) , X = Sl M2 2 , y = S \j2° ’ ( 6 ) 
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where p t denote the pion momenta^. To the same level of accuracy, the electromag¬ 
netic form factor is F(z) = 1 + z/ry with ry = MyjM\ ~ 2.5. In accordance with 
chiral counting, the polynomial in ([|) is assumed to have the general form 

Wf°\z)=a i + b i z (i — +,S) . (7) 


Up to a linear polynomial of this type, the dispersion integral corresponding to 
Fig. a is unambiguously calculable with the result 


wr(z) = — 


OLi + /?, 


Z~ Z 0 


F(z) x(z) , 


where z 0 = 1/3 + r 2 . The one-loop function 


4 4 r 2 i 4 r 2 

x(z) = ---^--(l-^)G(z/rl) 
9 3 z 


^4/z — 1 arcsin (y/z/2) 


G(z) = 


-\Jl-A/z In 




+ in 


1 + V 1 _ 4 /z 

satisfies y(0) = 0. In terms of the K —> 3n parameters in (|5j 
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a + = — ( 6 C + 62 ) , (3+ — 2 (d c + o^) > ots = —62 , 

The total form factor, given by 

Wi(z) = G F M 2 K ( ai + b iZ ) + Wr(z) , 


z <4 
z> 4 
we find 


fts — “ 3^2 • 


( 8 ) 


(9) 


( 10 ) 

( 11 ) 


is expected to be an excellent approximation to the complete form factor of 0(p 6 ). 
The polynomial piece has the most general form that can occur to this order, with a 
priori unknown low-energy constants contributing to the a+ 6 + The main assump¬ 
tion underlying the form factor (pT]) is that all other contributions to the dispersion 
integral except the two-pion intermediate state can be well approximated by a linear 
polynomial for small values of z. 

Using the central values of the K —*■ 37 t parameters |5| in Table [I], we can now 
calculate the branching ratios for K + —> n + £ + £~ and K$ —■► n°£ + £~ in terms of the 
corresponding parameters: 


B(K + - 

-> 7r + e + e' 

B(K + - 

> n + p + p 

B{K S - 

-> 7r°e + e 

B(K S - 

-> n°p + p 


0.14 - 3.23a+ - 0.886+ + 59.2a^ + 16.0a+6+ + 1.736 9 
1.13 - 19.2a+ - 6.326+ + 116a 2 + + 67.3a+6+ + 10.36*" 
0.01 - 0.76a s - 0.216s + 46.5a| + 12.9a 5 6 5 + 1.446|' 
0.07 - 4.52a s - 1.506 s + 98.7a| + 57.7a s b s + 8.956|" 


x 10" 8 
x 10~ 9 
x 10" 10 
x 10" 11 


J 


) 


( 12 ) 


3 The subscript 3 indicates the odd-charge pion, i.e. 7 r in K + —> 7r + 7r + 7r and 7r° in I\s 
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b c 

b 2 

d c 

d 2 

24.5 ± 0.3 

-3.9 ± 0.4 

-1.6 ± 0.3 

0.2 ± 0.5 


Table 1: Experimental values of K —> 37 r amplitudes || contributing to the IU/ r7r (z). 
All entries are in units of 10~ 8 . 


With a + and as expected to be of 0( 1), we observe that the amplitudes and 
branching ratios are actually dominated by the polynomial parts. W^ r ' K {z) in (J|) 
contributes to the rate mainly through its interference with (0). In fact, this really 
only applies to the K + mode because the two-pion contribution is practically neg¬ 
ligible for the neutral decay, due to the strong suppression of the K$ —■> 7 r + 7 r _ 7 r° 
amplitude. 


2.3 In the Brookhaven experiment BNL-E777 || both spectrum and branching 
ratio of the decay K + —> 7 r + e + e _ were measured. They fit a spectrum with two 
parameters C and A defined through 


dT 

dM ee 


C 


M ee M 


A 3 // 2 (l, ■ 


M 


M, 


K 


f 1 + A 



(13) 


where M ee is the invariant mass of the lepton pair (. M ee = M K y/z ), obtaining 


A = 0.105 ± 0.035 ± 0.015 , (14) 

B(K + -> vr+e+e") = (2.75 ± 0.23 ± 0.13) x 10 " 7 . (15) 

However, in order to subtract the background due to the process K + —> 7 r + 7 r°, 
7T° —> e + e _ 7, they make a cut in the dilepton invariant mass and consider only 
K + —> n + e + e~ events with M ee > 0.150 GeV. Thus the branching ratio actually 
measured is given by[] 


B(K + -»■ vr+e+e”) \ cut = (1.81 ± 0.17) x 10 “ 7 (16) 

and the result in (|1^) includes a theoretical extrapolation to the low M ee region. 

The branching ratio (|16|) can be translated into allowed domains for the param¬ 
eters a + and b + as shown in Fig. ^]. For a generous range \b + <2, two branches of 
solutions are in principle possible with opposite signs of a + . 

Unlike the rate, the spectrum is very sensitive to b + . The bounds on b + from the 
experimental spectrum are also indicated in Fig. ^]. Let us try to understand the 

4 We have extracted this value from ([L4]) and (|I^) using (|l3|) for the spectrum. For a conservative 
estimate of the error, we have scaled down the error given in (|l5|). 
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Figure 2: Allowed domains in the a + -b + plane (inside the dashed curves) from the 
experimental branching ratio ([H]). The bounds on b + from the measured slope (14) 
are also shown (full lines). 


allowed values of a + and b + with the help of the amplitude ([TT|). The real part of 
the two-pion contribution, Re W™(z), is a monotonically decreasing function of z 
in the physical region. Assuming first b + = 0, as would be the case at 0{jp 4 ) (except 
for the tiny contribution from the kaon loop), we infer that a + must be negative 
to reproduce the experimentally observed positive slope (0). In fact, the spectrum 
prefers a somewhat bigger slope requiring b + /a + > 0. This is made more explicit in 
Fig. |3;. where the theoretical spectrum (or rather the square of the form factor) is 
confronted with the experimental slope. 

At 0((p A ) only negative values of a + are allowed ||. However, in general there 
can be a second branch of solutions with positive a + . Since in this case the slope is 
negative for b + = 0, a sufficiently large b + (b + /a + > 1) is necessary to reproduce the 
observed spectrum. Analogously to Fig. |3], we exhibit the square of the form factor 
in Fig. []] for two sets of a+ and b + corresponding to the ±1<7 values of the measured 
slope. 

As we shall discuss in the next section, the second branch of solutions with 
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Figure 3: Comparison between different shapes of the Ii + —» 7 r + e + e~ form factor 
for a + < 0. The dotted and short-dashed lines correspond to the ±lcr values of the 
A parameter in (0. The solid and long-dashed lines correspond to the theoretical 
prediction for two sets of values of a + and b + . The comparison is to be applied for 
M ee > 0.150 GeV (vertical line) as measured by experiment 0. 


positive a + is disfavoured from a theoretical point of view. Indeed, naive chiral 
counting would suggest that the coefficient b + , which only arises at 0(p 6 ), should 
be smaller in magnitude than the leading coefficient a + . 

For the present analysis we do not discard the second type of solutions right 
away but confront the theoretical amplitude with another piece of experimental 
information that has recently become available through the measurement of the 
muonic decay mode of K + . The Brookhaven experiment of Adler et ah (7] has 
measured a branching ratio 

B(K + —*• 7r + / u + /i _ ) = (5.0 ± 0.4(stat) ± 0.7(syst) ± 0.6(th)) x 10 ~ 8 (17) 


which, together with the results from ||, implies a muon/electron ratio 


R = 


B(K + -> 7r + /i + /i') 
B(K + -»• 7 r + e + e") 


= 0.167 ±0.036 


(18) 
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Figure 4: Comparison between different shapes of the Ii + —> 7r + e + e form factor 
for a+ > 0 (as in Fig. |j). 

Can we use this experimental input to distinguish between the two branches in 
the a .|_ — b + plane? Looking first at the solutions with a + < 0, we find that R is a 
monotonically increasing function of |fq_| for all allowed values in the left branch of 
Fig. | and it is therefore always bigger than 0.23 (obtained for b + — 0). At first sight 
the second branch looks more promising for understanding the experimental result 
(p~8|). In fact, for b + = 0 the ratio R is small enough jljj to agree with flT8|). However, 
as Fig. |2] shows, the experimental spectrum for the electronic decay channel requires 
b + > a + . Since R is an increasing function of b + for a + > 0, the resulting values of 
R turn out to be even larger than in the previous case. Therefore, the theoretically 
unattractive solutions with positive a + cannot explain the small /i/e ratio either. 
For the values of a + , b + used in Figs. |3], f| the results are collected in Table [2|. 

On the basis of this model-independent analysis, we conclude that the central 
values of the /i/e ratio [0 and the slope of the electron spectrum |J together are not 
consistent. The ~ 2cr discrepancy can be due to a statistical fluctuation or could 
be an indication of peculiar non-standard physics. We stress that the inconsistency 
is essentially unrelated to the chiral expansion but only depends on the assumption 







a+ 

b + 

10 8 B (K+ ->• n + n+n~ ) 

R 

-0.68 

0.0 

6.78 

0.228 

-0.62 

-0.3 

7.29 

0.258 

0.55 

1.1 

7.19 

0.265 

0.47 

1.5 

7.89 

0.309 


Tabic 2: B(K + —>• n + p + p ) and the ratio R as a function of a + and b + . 


that both channels are dominated by the same K —> 7ry* form factor. Indeed, 
a form factor W + (z) rising with z as indicated by K + —> 7 r + e + e - data implies 
-R > -Rphase—space = 0.196 |]. 

3 Theoretical ideas on gl and b, 


The leading 0(p 4 ) predictions of a* and b, : are given by 


ai 4) = 


n (4) - 

CL q — 


^ ( 1/3 - ® + ) . 

Ltf 

-hi(i/ 3 _» s ), 

Ltf 


bf - 

J4) = 


G f 60 ’ 

Gg J_ 

G f 30 ’ 


(19) 


in terms of the usual parameters Gg and Wi [|l|. Local contributions to the bi are 
forbidden at 0(p 4 ) and the small values in ([T9]) are generated by the expansion of 
the kaon-loop function. Sizable corrections to the bi are expected at 0(p 6 ) where 
local terms are allowed. However, since local terms contribute to the a* already 
at 0(p 4 ) (Wi terms), naive chiral counting would suggest ~ 0(p 6 ) fO(p 4 ) < 1. 
The solution with both a + and b + negative discussed in the previous section satisfies 
this expectation. 

The expressions of w + and ws in terms of the 0(p 4 ) low-energy couplings 


10], Lg [11] are 


w + = 


Ws = 


647r 2 i / u 2 

— [NIM -Nl M + ZF M] + -In 

{2NM + AT[ 5 M] + lln 1 I 


Ml, 


( 20 ) 
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Unfortunately, the values of the N[ are not known and to make definite predictions 
we need to rely on model-dependent assumptions. For this purpose, it is interesting 
to note that while the combinations of low-energy constants appearing in w + and Wg 
(or a + and as) depend separately on the renormalization scale /i, the combination 
occurring in a + + as does not. As originally noted in |I|], this scale independence 
might be related to the structure of the effective four-fermion Hamiltonian relevant 
for K —> ni + £~. Indeed, as we will discuss in the next section, there is only one 
dimension-six operator giving a non-vanishing contribution to the a* at leading 
order: Q^y = S7 M (1 — 7 , 5 ) 0 ? Ay,A- Due to the octet structure of the s^d current, this 
operator does not affect the combination a + + ag- Actually, the cancellation of this 
short-distance contribution holds not only for a + + as but (in the limit of isospin 
conservation) for A(K + —> n + £ + i~) + A(Kg —> tt°£ + £~) in general, that would thus 
be completely determined by low-energy dynamics. As a consequence, we find it 
more reliable to predict the sum rather than the separate expressions of W+° l and 
Wg° l using a low-energy model. In the following we shall employ the Vector Meson 
Dominance (VMD) hypothesis to estimate W + o1 + W| o1 . 

The vector meson contributions to the low-energy constants of the 0(p 4 ) weak 
Lagrangian have been discussed in ||T 0 | , |T 2 | . Employing the vector field formulation 
of Ref. [ 121, the result for Wy° x -{-Wg 01 is independent of the factorization hypothesis 
and can be specified in terms of a single unknown parameter rjy 


wH) _L_ w(4) _ 

VV +,V T Wgy — 

Lrp L 


16vr 2 fy ( 2r) V - 1 ) + - In 


M, 


M: 


K 


( 21 ) 


where \ fy\ ~ 0 . 20 , as obtained from T(p° 


e T e' 


"). The logarithmic term in 


is 


a residual effect of the loop amplitudes, whereas the term proportional to fy is the 
local vector meson contribution. The rjy parameter is not known but is expected to 
be in the range 0 £ rjy < 1 ; in principle it is measurable in other processes [0- Note 
that for typical values of r/y (but for rjy ~ 1 / 2 ) the local vector meson contribution is 
dominant. A relation equivalent to © was obtained in Jl(J under the factorization 
assumption (the result of [[TTJ] is obtained from (|2T|) with the substitution rjy —> kp)- 
j! 0 [], the separate predictions of as and a + involve an additional 


As discussed in 


unknown coupling which is difficult to interpret in the framework of VMD. 

The evaluation of the 0(p 6 ) local terms generated by vector meson exchange is 
more involved [|13|, |0. In general, we can identify two kinds of contributions: those 
originating from genuine 0(p 6 ) weak transitions and those generated by the pole 
expansion of the leading 0(p 4 ) VMD results (we refer to O, [Ll for a detailed dis¬ 


cussion). Interestingly enough, all the effects of the genuine 0(p 6 ) weak transitions, 
evaluated under the factorization assumption, drop out in the sum W+° l + Wg° l . 
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Moreover, the 0{p 6 ) vector exchange contributions only modify the slope parame¬ 
ters bi, leading to the simple result 

Gs i. 


W® + W)?> = 


-y 


( 6 ) _ 


sy 


G f r 'y L 


167T ./^(2 W - 1) 


( 22 ) 


Combining (|2l|) and (|22|) and neglecting the presumably small logarithmic term in 
(0. we then obtain 

X - ( 23 ) 


H 6 ) 4 - //' 

°+y + °sy 


l(6) 


,(4) 


,(4) 


, y + a s'v 1 v ' 

In order to disentangle neutral and charged channels beyond the above relation 
one has to rely on additional assumptions. We notice that out of the two solutions 
for a + and b + discussed in the previous section, the one with a + < 0, consistent 
with chiral counting, is also compatible with the relation b + /a + = 1 jry. This would 
follow directly from the assumption of vector meson dominance in the polynomial 
part of the form factor, i.e. from the hypothesis 


Wr\z) = a 


V 


Ty ~ * 


(24) 


As commented previously, both form factors Wy° l and Hd? 01 may receive short- 
distance type contributions that in principle could spoil this relation. However, 
the form factor of the matrix element (7r|s7 M d| K) exhibits a clear K*( 892) pole 
dominance (in the SU(3) limit it coincides with the electromagnetic form factor 
F(z) discussed in the previous section). Hence, we conclude that the pole structure 
of Wl° ] (z) is not restricted to the long-distance part. Only the genuine 0(p 6 ) weak 
vector transitions could spoil the relation (^4|). If these were negligible, and in 
general this is not necessarily the case, the ratio b^/di — 1/r'y would hold for both 
channels. 

We can use this assumption to analyze the Kg —■► 7r°£ + £ _ mode. Taking a look 
at the expressions of B(Ks —^ n°£ + £~) in (|12j) and considering that asbs > 0 (as 
follows from bg/as = 1 /r'y), we observe that for \as\ ~ 0.2 (i.e. for B(Kg —^ 
7r°e + e _ ) > 2 x 10“ 10 ) the constant and linear terms in ag and bg are negligible and 
can be dropped. In this way we obtain an expression for the branching ratios in 
terms of just one parameter 


B(K S —■> 7r°e + e ) ~ 5.2 a 2 s x 10 9 

T3(TS - ^ 1 O ^ ^ 1 n—9 


B(K S —7T° jji + p ) ~ 1.2 a 2 s x 10 


from where we predict 


B(Kg 7TV+/X-) 
B{K S -»■ 7r°e + e-) 


~ 0.23 . 


(25) 


(26) 
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Unfortunately the lack of information on ag does not allow us to predict the separate 
rates. We note however that if as ~ 0(1), as expected on general grounds, the 
electron mode should be within reach of the KLOE experiment 


4 CP violation in K —> n£ + £ decays 


The dominance of single photon exchange does not apply to the UP-violating parts 
of K —> n£ + £~ amplitudes. In this case the hierarchy of the CKM matrix implies 
that Z° -penguin and fU^-box diagrams, dominated by short-distance contribu¬ 
tions, play an important role. As a result, UP-violating amplitudes can be more 
conveniently studied by means of an appropriate four-fermion Hamiltonian that, in 
the case of s —> d£ + £~ transitions, is known to next-to-leading order ]15|] . At scales 
H < m c it is given by 


H 


|AS|=1 

eff 


G f 


6,7V 


y/2 


Ks v ud[ J2 CO + r yiO)QiO + ry 7A (M w )Q 7A {M w ) 


i=1 


+ h.c., 
(27) 

where r = — (V t *V t d) / (V* s Vud) and denote the CKM matrix elements. Here Q i ;2 
are the current-current operators, Q 3 ,..., Q 6 the QCD penguin operators, whereas 

Qiv = S 7 M (1 - and Q 7A = s^(l - (28) 


are generated by electroweak penguin and box diagrams (0 |T7| . In the standard 
phase convention, the overall factor V* s V u d is chosen to be real and direct CP viola¬ 
tion is driven by the imaginary part of r. With this choice of phase, only the QCD 
penguin operators and Q 7 y , Q 7A are relevant for estimating direct UP-violating am¬ 
plitudes. As can be understood by looking at the corresponding matrix elements, 
the dominant role is played by Q 7 v,Qia HH, with Viv (1 GeV) ~ 5.7 x 10~ 3 and 
y 7A (M w ) ~ —5.3 x 1CT 3 JTH) (corresponding to m t {m t ) = 167GeV). 

The contribution of Q 7 y to K decays interferes with the long-distance 

amplitude discussed in the previous sections. In fact, the contribution proportional 
to the coefficient z 7 y is already included in the polynomial part of the K —> 777 * 
form factor (as discussed in section |3|). On the contrary, the interference of Q 7A 
with the photon exchange amplitude vanishes as long as the lepton polarizations are 
summed over. In the following we shall discuss the role of y 7 vQ 7 V hr Kl 
and K ± —> vUe+e - decays. 


7r°e + e 


4.1 From the definition ([[]), UP invariance would imply that the K L —> 7r°7* form 
factor Wl{z) vanishes. In the limit of UP conservation, the decay Kl —> 7r°e + e _ can 
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7T°77 


only proceed through the two-photon process K L 
or via subleading terms in the expansion of the W ± -box diagram 
mechanism is expected to be dominant yielding B(Kj, 


7 r°e + e 


19 


21 


The former 
7r°e + e - )cpc ~ few x 10' 12 


13, 20 


In the presence of CP violation, the Kl —> 7r°e + e _ transition can also proceed 
through the small £kK\ piece of the K L wave function (indirect CP violation) or via 
the direct (TP-violating amplitude generated by (p7|) . The contribution of yjyQiv 
appears in the function IT]) 01 , defined in analogy to (0), with 


47r Im X t 1 

a L =-i—j= - Dtv and — — ~r •> 

V2 a cll r v 


(29) 


where A t = VtdV* s . On the other side, the indirect C’P-violating contribution is given 
in terms of the parameter as A(K L —> 7r° e + e ) ind — \eK\e llT A(K s —> 7r°e + e ). 
Summing the two vector-like OP-violating contributions, one has 


wr 1 = 6.9 x 10 


-4 


3.3 as e* 4 - i- 


. Im X t 


10 


-4 



2 


1 H- 2~ 


1 r v\ 


(30) 


Unfortunately the sign of as is not known. Moreover, we recall that for |cis| £ 0.2 
the dominant term in B(Kg —> 7T°e + e _ ) is provided by a|, therefore it is hopeless 
trying to extract information about the sign of as from the measurement of B(Ks 
7 r°e + e _ ).[] However, since as is real a strong interference (constructive or destructive) 
is expected for |as| and |Im A*/10 —4 1 of the same order of magnitude. 

Including the C’P-violating contribution of the Qta operator that does not in¬ 
terfere with ©, we collect all the C’P-violating terms in Kj j —> 7i 0 e + e~ obtaining 


B(K l 


7r°e H 


)CPV 


15.3 a 2 s 


f. o Im A t , „ $ 

6.8 4 a s + 24 


Im A t 

i n- 4 


X 10 


-12 


(31) 


From (^) and (|3lD a very interesting scenario emerges for as ~ —0.5 or as ~ 1.0. 
Since Im X t is expected to be ~ 10” 4 , one would have B(K L —»■ 7r°e + e _ )cpy ~ 10~ n 
in this case. Then the CP-conserving contribution, which does not interfere with 
the C'P-violating part in the rate, could be neglected. Moreover, the Ks —► 7r°e + e” 
branching ratio would be large enough to allow a direct determination of \as\- Thus, 
from the interference term in ([U]) one could perform an independent measurement 
of Im X t , with a precision increasing with the value of \as\- 

5 In principle the relative sign of as with respect to A{Kg —* n + Tt~C) could be measured by 
a careful analysis of the Ks —> 7r°e + e _ spectrum (similarly to a+ in the charged mode). However, 
this program is made very difficult by the smallness of A(Ks —> 7r + 7r _ 7r°). In addition, even if it 
were possible, one would then have to rely on a model-dependent assumption about the sign of 
A(Ks —^ 7r + 7r _ 7T°) from (p7|). 
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A relation similar to (31) can be traced from the work of Dib et al. [|17| (see 
also |TE|). We note that this result arises from the VMD assumption about the 
K s —■► 7r°e + e _ form factor discussed in section [| which leads to (^). According to 
this hypothesis, the relative weight of direct and indirect C'P-violating components 
in ([H]) is not affected by possible cuts in the z spectrum. On the contrary, we 
recall that the CP-conserving rate does have a different (softer) spectrum in z: 
essentially the factor A 3 / 2 (l, z, r 2 ) in @ is replaced by A 5 / 2 (l, z, r 2 ) Jl9], |2l[, |22[]. This 


fact could provide an additional handle for disentangling the various components of 
K l 


7r°e + e 


4.2 The interference of the long-distance K —> 7T7* amplitude and the short- 
distance contribution of yzvQzv leads to an asymmetry between the widths of K + —> 
7r + e + e~ and K~ —> 7r _ e + e _ , which is a clear signal of direct CP violation [|T9|. This 
quantity, defined by 


T(A' + -»• 7r + e + e 
/■(l —Or ) 2 


T/T'- - + -\ GpOl M k 

— T A —> 7 r e e ) = -f— x 

V ; 768 vr 5 

2 \ 


/4 r 2 


dz A 3/2 (l, r 2 ) y 1 - 4^ I 1 + 2^ J Im Wl°\z) Im W™(z) , 


can be calculated unambiguously up to a sign. Indeed, the only contribution to 
Im Wl° l comes from the short-distance regime and it is given by 


47r Im A* 

Im a + = —= - y 7V 

V2 a 


Im b + 
Im a + 


„2 ’ 
v 


(32) 


whereas the imaginary part of W™ is determined by the physical K + —> 7r + 7T + 7r _ 
amplitude. Using the results in (|8- 1C) and the experimental value of T(A' + — 
7T + e + e~), we get 


|T(A' + —> 7r + e + e ) — Y{K n e + e 


T(K + —»■ 7r + e + e ) + T(AT —> 7r e + e ) 


= (0.071 ±0.007) x |Im A t | . (33) 


The numerical value in (|33|) corresponds to an increase of a factor ~ 2 over the 
estimate at leading order in the chiral expansion [ TO |. However, given that Im X t ~ 
10 -4 , it is still very difficult to detect such an effect, at least within the Standard 
Model. 

A more interesting observable is the unintegrated asymmetry defined by 


ST(z) = 


— (A' + —> 7r + e + e ) — —{K —»■ n e + e ) 
dz dz 


T(A' + —> 7r + e + e ) ± T(K —> n e + e ) 


(34) 
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Figure 5: Differential charge asymmetry defined in (]34|). The dashed line is the 
leading-order result in the chiral expansion, whereas the full line includes the uni- 
tarity corrections. 


that we have plotted in Fig. both at leading order in the chiral expansion and 
including unitarity corrections. This asymmetry is seen to have a maximum around 
z ~ 0.4 that, being far away from the dominant background, could be more easily 
accessible from the experimental point of view than the integrated asymmetry. 


5 Conclusions 

We have performed a model-independent analysis of K + —> 7t + DT - and Kg —■> 
7 t°£ + K decays including K —> 37T unitarity corrections and a general polynomial 
decomposition of the remaining dispersive amplitude. Using as input the presently 
available data on Ii + —> 7r + e + e _ [jSJ, we concluded that the ratio R = B(K + —> 
7 t + / j, + B(K + 7r + e + e ~) must be larger than 0.23 unless there is a contribution 
from some non-standard physics. This result is to be compared with the recent 
measurement R = 0.167 ± 0.036 0. 

We have shown that it is very difficult to predict B(Kg —> 7r°e + e _ ) with- 


15 




out strong model-dependent assumptions, except for an approximate upper bound 
B(Ks vr°e + e _ ) < 10~ 8 . On the other hand, a realistic estimate of the muon/elec¬ 
tron ratio can be obtained in the framework of VMD: B(Ks —■► n°/u, + fi~)/B(Ks —■► 
7T°e + e~ ) ~ 0.23. 

We have reanalysed the CP-violating contributions to Kl —> 7r°e + e _ , empha¬ 
sizing a possible strong interference (destructive or constructive) between direct 
and indirect CP violating amplitudes. Together with a precise determination of 
B ( Ks —► 7r°e + e _ ), this could provide a possible handle for isolating the direct CP- 
violating component. Finally, while the unitarity corrections lead to an increase 
of about 100% for the charge asymmetry T(K + ii + e + e~ ) - T(K~ —> 7r _ e + e _ ), 
the Standard Model prediction is still too small to be within reach of forthcoming 
experiments. 


Note added 

At the recent ICHEP98 in Vancouver, preliminary results from a high-statistics 
K + —> 7r + e + e~ experiment (BNL-E865, presented by Hong Ma) were announced. 
While the mean value of the decay rate agrees with the previous measurement ||, 
the slope is significantly bigger: A = 0.20 ± 0.02. Although such a big slope is still 
compatible with the theoretical expectation |6 + | < |a + |, the parameter |6 + | would 
have to be bigger than the naive VMD prediction b + = a + /ry , but consistent with 
the factorization model predictions for the genuine 0{p 6 ) vector meson contribution. 
Finally, we observe that the larger slope aggravates the discrepancy between theory 
and experiment for the /i/e ratio R. 
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